Dynamics of Atom-Mediated Photon-Photon Scattering I: Theory 
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The mediated photon-photon interaction due to the resonant Kerr nonlinearity in an inhomo- 
geneously broadened atomic vapor is considered. The time-scale for photon-photon scattering is 
computed and found to be determined by the inhomogeneous broadening and the magnitude of the 
momentum transfer. This time can be shorter than the atomic relaxation time. Effects of atom 
statistics are included and the special case of small-angle scattering is considered. In the latter case 
the time-scale of the nonlinear response remains fast, even though the linear response slows as the 
inverse of the momentum transfer. 



I. INTRODUCTION 

Recently there has been experimental and theoretical 
interest in the nonlinear optics of confined light A 
medium possessing an optical Kerr nonlinearity and con- 
fined within a planar or cylindrical Fabry-Perot resonator 
gives rise to new nonlinear optical phenomena such as 
soliton filtering and bilateral symmetry breaking pLp|. 
The classical nonlinear optics of this system is described 
by the Complex Ginzburg-Landau equation (CGLE) 
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where E is the electric field envelope, k is the longitudinal 
wavenumber, u> — ck/riQ is the field envelope angular fre- 
quency, A is a mode overlap factor, Afc is the wavenum- 
ber mismatch from the linear-cavity response and T is 
the field amplitude decay rate. The classical dynamics 
of Eq. (0) describes the mean-field behavior of a system 
of interacting photons coherently coupled to an external 
reservoir. A photonic system of this sort is a versatile 
model system for condensed matter physics in reduced 
dimensions [Q, as the parameters Ak,ri2, I\ and E^ in 
Eq. (|l|) are subject to experimental control. In partic- 
ular, an atomic vapor can provide a strong Kerr nonlin- 
earity which is tunable both in strength and in sign. In 
this case the nonlinearity arises from the saturation of 
the linear refractive index, which is a strong function of 
the drive laser frequency near an absorption resonance. 

Some of the most interesting proposed experiments 
for this system, including generation of few-photon 
bound states || , direct observation of the the Kosterlitz- 
Thouless transition in an optical system Q] and observa- 
tion of quantum corrections to the elementary excitation 
spectrum of a ID photon gas |j]|7]] intrinsically involve 
photon correlations. For this reason, it is important to 
understand the microscopic (and not just mean-field) be- 
havior of photons in an optical Kerr medium. We specifi- 
cally consider saturation of the resonant electronic polar- 
ization of a Doppler-broadened atomic vapor, a medium 



which has been proposed for quantum cavity nonlinear 
optics experiments and used to observe a nonlinear cavity 
mode |Q . Thus the system under consideration involves 
dispersion, loss, inhomogeneous broadening, and the con- 
tinuum of transverse modes in an extended resonator. 

Sophisticated techniques have been developed for 
treating mediated interactions among photons in non- 
linear media. One approach is to obtain an effective 
theory in which the quanta are excitations of coupled 
radiation-matter modes, by canonical quantization of the 



or by direct attack 
This approach has 



macroscopic field equations 
on a microscopic Hamiltonian 
the advantage of generality and is suited to multi-mode 
problems, but has basic difficulties with loss and disper- 
sion near resonance pd|-p^[|. Microscopic treatments in- 
clude Scully-Lamb type theory [^4|Jl5| and application of 
phase-space methods [ p|[L7| |. A strength of these tech- 
niques is their ability to handle relaxation and popula- 
tion changes. They are, however, cumbersome to apply 
to inhomogeneously broadened media and to multi-mode 
problems. 

In this paper we characterize the atom-mediated 
photon-photon interaction using an accurate microscopic 
model and perturbation calculations. This allows us to 
determine the time-scale of the mediated photon-photon 
interaction in the atomic vapor, despite the complexity 
of the medium. We find that the interaction is fast and 
not intrinsically lossy, even for small momentum trans- 
fer. Thus the medium is suitable for quantum optical 
experiments, including experiments using the NLFP as a 
model for the interacting Bose gas. 



II. SCATTERING CALCULATIONS 

The complete system is treated as the quantized elec- 
tromagnetic field interacting via the dipole interaction 
with an vapor of atoms of mass M. The perturbation 
calculations are performed in momentum space, as is nat- 
ural for thermodynamic description of the atomic vapor. 
This also makes simple the inclusion of atomic recoil ef- 
fects. The dipole interaction term is identified as the 
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perturbation, so that the eigenstates of the unperturbed 
Hamiltonian are direct products of Fock states for each 
field. In the rotating wave approximation, the unper- 
turbed and perturbation Hamiltonians are 
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where ak, Q is the annihilation operator for a photon of 
momentum hk. and polarization a, c„ iP is the annihila- 
tion operator for an atom in internal state n with center- 
of-mass momentum hp, E is the quantized electric field 
and d is the atomic dipole field. Polarization plays only a 
very minor role in this discussion so polarization indices 
will be omitted from this point forward. 

The simplest mediated interaction is photon-photon 
scattering, which transfers momentum from one pho- 
ton to another by temporarily depositing this momen- 
tum in the medium. Specifically, photons with mo- 
menta k, 1 are consumed and photons with momenta 
k' = k + q, 1' = 1 — q are produced. The lowest-order 
processes to do this are fourth order, so we look for rel- 
evant terms in H' H 1 H' H 1 . A parametric process, i.e., 
one which leaves the medium unchanged, sums coher- 
ently over all atoms which could be involved 18 1. Due 
to this coherence, the rates of parametric processes scale 
as the square of N/V, the number density of atoms. In 
contrast, incoherent loss processes such as Rayleigh and 
Raman scattering scale as N/V. Thus for large atomic 
densities, a given photon is more likely to interact with 
another photon than it is to be lost from the system. 

In this sense, the interaction is not intrinsically lossy, 
as are some optical Kerr nonlinearities such as optical 
pumping or thermal blooming. The latter processes re- 
quire absorption of photons before there is any effect on 
other photons. For this reason, they are unsuitable for 
quantum optical experiments such as creation of a two- 
photon bound state. 

One parametric process, photon-photon scattering at 
a single atom, is described by the diagram of Fig. 1. The 
relevant terms in H'H'H'H' contain 



C Lp C rf,P+l' a l' C d,p+l' C c,p-q a l 
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or permutations k' <-» 1', k <-> 1 for a total of four terms. 
Here p is the initial atomic momentum and a through d 
index the atomic states involved. With the assumption 
that no atoms are initially found in the upper states b 
and d, i.e., n b = n d = 0, this reduces to 
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where the n are number operators for the atomic modes 
and the upper and lower signs hold for Bose and Fermi 
gases, respectively. The difference for atoms of differ- 
ent statistics reflects the fact that the scattering process 
takes the atom through an intermediate momentum state 
which could be occupied. Occupation of this intermedi- 
ate state enhances the process for Bose gases but sup- 
presses it for Fermi gases. 

A thermal average of the relevant terms in H'H'H'H' 
gives the thermally averaged effective perturbation 
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where 
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and similar expressions obtain for . and 

the average occupancy of the atomic state \a,p). The 
R; are the resonance denominators 
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Here hu>ij = h(tOi — u>j) is the energy difference between 
states i and j, 7, is the inverse lifetime of state i and rj 
is a vanishing positive quantity. Here and throughout, 
the process is understood to conserve photon momen- 
tum, but for clarity of presentation this is not explicitly 
indicated. 

As described in Appendix [A|, intensity correlation func- 
tions for photon-photon scattering products contain a 
Fourier transform of the scattering amplitudes 
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where 6k' is the output photon energy shift, xa,b and 
iA.B are detection positions and times, respectively, and 
t_ = £b — %b/c — tA + ^a/c is the difference in retarded 
times. This expression allows us to determine the time 
correlations for photon-photon scattering in a number of 
important cases. 



III. LARGE-ANGLE SCATTERING 

The simplest configuration to understand is that of 
counterpropagating input beams producing counterprop- 
agating output photons scattered at large angles. This is 
also the most convenient experimental geometry. 



A. One Atom Process 

Scattering amplitudes and rates for right-angle scatter- 
ing by the one-atom process are shown in Fig. 2 and Fig. 
3, respectively. For the moment we ignore the statistical 
correction due to the « Q ,p"e,p-q term in Eq. (^), which 
will be considered separately. The the vapor is treated as 
a gas of two-level atoms. The parameters are the Doppler 
width (5d = M&bT/M) 1 / 2 , where &b is Boltzmann's con- 
stant, the radiative linewidth 7^ = As/2 where Ae is the 
Einstein A coefficient, and the detuning A = ck — ujba, in 
the ratios 7b = 0.01<5d, A = 2tt6b- The amplitude units 
are arbitrary, but do not vary between graphs. 

At this point it is important to note that the dura- 
tion of the correlation signal is much shorter than the 
coherence lifetime of an individual atom, approximately 
% x . In fact, the duration of the correlation signal is de- 
termined by the momentum distribution, a property of 
the medium as a whole. This can be explained in terms 
of the coherent summation of amplitudes for scattering 
processes occurring at different atoms. The process is 
coherent only when it is not possible, even in principle, 
to tell which atom participated. This clearly requires 
momentum conservation among the photons, but it also 
limits the duration of the atomic involvement. An atom 
acting as intermediary to transfer momentum q is dis- 
placed during the time in remains in the state c of Fig. 1. 
If this displacement is larger than the thermal deBroglie 
wavelength A it is possible, in principle, to determine 
which atom participated. This limits the duration of the 
coherent process to 8r ~ AM/frq. 



B. Statistical Correction 

As noted above, the quantum statistics of the atoms 
in the vapor contribute a correction to the single-atom 
scattering amplitude. This correction (with the sign 
appropriate for Bose atoms) is shown in Fig. 4 for 
a gas with phase space density NA 3 /V = 1/2, where 



A = (Mk B T/2TTh 2 ) 1 / 2 is the thermal deBroglie wave- 
length. Parameters are as for Fig. 2. 

C. Simultaneous scattering 

A second parametric process, simultaneous scattering, 
is described by the diagram of Fig. 5. The relevant terms 
in H'H'H'H' contain 

4,p c d,P+i' a [' c I, P -q c f),P+kaL 

XC d,p+l' C c,p-q a l C b,p+k C a,P a k (12) 

or permutations k' *-* V, k <-» 1 for a total of four terms. 
Making the same assumption as before, this reduces to 

^a.p^cp-qa^a^aiak- (13) 

This process corresponds to the absorption of each pho- 
ton by an atom before emission of either, and thus de- 
scribes a two-atom process and is of the same order in 
the atomic number density as the Fermi and Bose cor- 
rections to single-atom scattering. The kinematical and 
geometric factors of Eq. (Q) and Eq. (||) are the same for 
this process, and the resonance denominators are 

flf) = c (k + I - k') - A[p . 1' + Z'2/2] _ uj da + i ld 

4 2) = c(k + l) - A[p . k + k 2 /2 + (p - q) • 1 + Z 2 /2] 

-UJ ba + ijb - Lddc + ijd 

4 2) = c(fc) - -^[p • k + k 2 /2] - u ba + i lb . (14) 

Amplitudes for simultaneous scattering are shown in 
Fig. 6 for a gas with a phase space density of one half. 
Parameters are as for Fig. 2. 

D. Fermi and Bose Gases 

The statistical correction and two-atom scattering con- 
tributions add coherently, giving considerably different 
correlation functions for moderate degeneracy Bose vs. 
Fermi gases. This is illustrated in Fig. 7 and Fig. 8, 
which show the scattering rates vs. delay for Bose and 
Fermi gases with a phase space density of one half. Pa- 
rameters are as for Fig. 2. 

E. Ladder Process 

In atoms with a "ladder" level structure, in which three 
levels a-c are ordered in energy ui c > u b > u> a and con- 
nected by matrix elements /i& a ,Mcb 7^ 0, fi ca = 0, an 
additional process described by the diagram of Fig. 9 is 
possible. The relevant terms in H'H'H'H' contain 
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c i,p c rf,P+i' a i' c d,p+i' c c,P+k+ia k ' 

X c lp+k+l c fc,P+kaici,p+k c a,P a k (15) 

or permutations k' <-> 1', k «-> 1 for a total of four terms. 
Making the same assumption as before, this reduces to 

n a,p a i' a k' a i a k- (16) 

This process corresponds to the absorption of both pho- 
tons by an atom before emission of either, and thus de- 
scribes a one-atom process which is of the same order in 
the atomic number density as one-atom scattering. The 
kinematical and geometric factors of Eq. (fy Eq. ([)[) are 
the same for this process, and the resonance denomina- 
tors are 

R{ 3) = c(k + I - k') - A[p . i' + J'2/2] - uj da + i ld 

4 3) = C (fc + 0-|:[p-(k + l) + |k + l| 2 /2] 

-uj ca + ij c 

4 3) = c(k) - A [p . k + k 2 /2] - u ba + i lh . (17) 
Ivl 

Right-angle scattering amplitudes for this process are 
shown in Fig. 10. Parameters are as for Fig. 2. 

F. Lorentz-model Behavior 

It is interesting to consider the case of a ladder atom 
with equal energy spacing lu c i, — u>b a and matrix elements 
\fj. c b\ 2 = 2|^i;,a| 2 - In this case the states a-c are equivalent 
to the lowest three levels of a harmonic oscillator, i.e., to 
a Lorentz model, and the medium is effectively linear for 
two-photon processes. 

The amplitudes for the one atom process of Eq. (Q) 
and the ladder process of Eq. (|l5|) partially cancel. The 
resulting signal is smaller and lacks oscillations, as shown 
in Fig. 11. Parameters are as for Fig. 2. 

G. Background Events 

In addition to the photon-photon scattering processes, 
Rayleigh scattering (and Raman scattering for more 
complicated atoms) will create an uncorrelated coinci- 
dence background. This background is calculated in Ap- 
pendix |A|. The coincidence signal, consisting of both 
the Lorentz-model atom photon-photon scattering sig- 
nal and the incoherent background is shown in Fig. 12. 
The peak coincidence rate (at St = 0) is approximately 
twice the background, accidental coincidence rate. In 
the limit of large detuning, it becomes exactly twice ac- 
cidental rate. This can be explained in analogy with 
the Hanbury-Brown-Twiss effect as follows: For the op- 
timal geometry the drive beams are conjugates of each 
other H(x) = G*(x) and the detectors are in opposite 



directions. The linear atoms act to create a random in- 
dex grating which scatters a chaotic but equal (up to 
phase conjugation) field to each detector. As expected for 
chaotic light jl!| , the fourth-order equal-time correlation 
function is twice the product of second-order correlation 
functions. 

(£ 2 (x A ,t)£ 2 (x B ,t)) = 2(E 2 (x A ,t)) (E 2 {* B ,t)) . (18) 



IV. SMALL-ANGLE SCATTERING 

Thus far the discussion has involved only large-angle 
scattering. In the context of cavity nonlinear optics all 
fields are propagating nearly along the optical axis of the 
cavity so it is necessary to consider scattering processes 
for nearly co-propagating or nearly counter-propagating 
photons. As argued above, the temporal width of the 
correlation signal scales as 1/g, the inverse of the mo- 
mentum transfer. This is shown in Fig. 13 and Fig. 14, 
which show rates for scattering photons from beams in 
the x—z plane into the the y—z plane. In all cases the 
beam directions are 0.1 radian from the z axis. The co- 
incidence distribution shows oscillations which die out 
on the time-scale of the inverse Doppler width, and a 
non-oscillating pedestal with a width determined by the 
momentum transfer q. 

The pedestal, however, does not correspond to the du- 
ration of the nonlinear process in this case. As above, by 
considering a ladder atom with the energy spacings and 
matrix elements of a harmonic oscillator we can isolate 
the linear optical behavior. As shown in Fig. 15 and 
Fig. 16, this behavior includes the pedestal, but not the 
oscillations, indicating that the nonlinear optical process 
is still fast, with a time-scale on the order of the inverse 
Doppler width. 



V. LIMITATIONS ON SCATTERING ANGLE 

Due to the limited width of the atomic momentum 
distribution, the resonance denominator is small if 
the input and output photons are not of nearly the same 
energy. Since the complete process must conserve pho- 
ton momentum, input photons with net transverse mo- 
mentum in the output photon direction will scatter less 
strongly. The width of this resonance is very narrow: a 
net transverse momentum k y + l y ~ k\/ k^T/Mc 2 is suf- 
ficient that few atoms will be resonant. As ^Jk-gT /Mc 2 
is typically of order 10~ 6 in an atomic vapor, this would 
be a severe restriction on the transverse momentum con- 
tent of the beams in a cavity nonlinear optics experiment. 
However, as shown in Fig. 16, the narrow resonance as- 
sociated with contributes the linear response of the 
medium. The nonlinear response, which has the same 
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resonance character as the "ladder" process, is not lim- 

(3) 

ited in this way because R 2 does not depend upon the 
output photon energies. 



VI. OUTPUT POLARIZATION 

The polarization of the output photons depends on 
the structure of the atom and can produce polarization- 
entangled photons. For example, if the input photons are 
propagating in the ±z directions and are x polarized, the 
two absorption events in the above diagram change the z 
component of angular momentum by 5m = ±1. In order 
for the process to return the atom to its initial state, the 
two emission events must both produce 5m = ±1 or both 
5m = 0. For right angle scattering with the detectors in 
the ±y directions, the output photons must therefore be 
either both x or both z polarized. If both polarizations 
are possible, the emitted photons are entangled in polar- 
ization, as well as in energy and in momentum. 



VII. CONCLUSION 

Time correlations in photon-photon scattering provide 
an indication of the time-scale over which the atomic 
medium is involved in the interaction among photons 
in a nonlinear medium. It is found that the time-scale 
is determined by the inhomogeneous broadening of the 
medium and the magnitude of the momentum transfer. 
For large-angle scattering, the time-scale of involvement 
is 5t ~ AM/hq, while for small-angle scattering the time- 
scale is 5t ~ AM/hk. As this time-scale is shorter than 
the atomic relaxation time, calculations which contain an 
adiabatic elimination of the atomic degrees of freedom 
necessarily overlook the fastest dynamics in this process. 



APPENDIX A: PHOTON CORRELATIONS 

1. Detection Amplitudes 

Unlike a genuine two-body collision process, atom- 
mediated photon-photon scattering has a preferred ref- 
erence frame which is determined by the atomic mo- 
mentum distribution. To calculate the photon cor- 
relations we work in the "laboratory" frame and as- 
sume the momentum distribution is symmetric about 
zero. We consider scattering from two input beams 
with beam shapes G(x) = y -1 / 2 Ek^C 4 ) exp[ik • x] and 
if (x) = V~ x l 2 h(l) expfil • x] which are normalized as 
Ek \9( k )\ 2 = Ei \ h 0)\ 2 = !■ We further assume that the 
beams are derived from the same monochromatic source 
and are paraxial, i.e., that <?(k) is only appreciable in 
some small neighborhood of the average beam direction 
ko, and similarly for h(V) around L> The geometry is 
shown schematically in Fig. 17. For convenience, the 



beams are assumed to each contain one photon, so that 
the initial state of the field is 



I0(o)) = 44 |o> 



(Al) 



where the creation operators 4;>4r are = 
^ k <7(k)aj c and A^ H = Ei^(^) a r Scaling of the result 
to multiple photons is obvious. 

We use Glauber photodetection theory to determine 
the rates at which scattering products arrive at two detec- 
tors A and B at space-time points (x^, t A ) and (x#, <#), 
respectively. We compute the correlation function in the 
Heisenberg representation 



P(xA,tA,X B ,t B ) 



4 +) (x B ,t B )4 +) (xA,tA)|0(O)) H | 2 

where the photon field operator is 

$W (x, t) = V- 1 ' 2 J2 a Ka (t) exp[ik • x]. 



(A2) 



(A3) 



k . a 



This field operator is similar to the positive fre- 
quency part of the electric field and is chosen so that 
$( - )(x, t)<E>( + )(x, t) is Mandel's photon-density operator 
p5] |. To make use of perturbation theory Eq. ([A|) is 
more conveniently expressed in interaction representation 
as 

P(x A ,t A ,X B ,t B ) 

= | <0| $i +) (x B ,i B )[/ / (i s ,iA)$| +) (xA,iA) \4>{t A )) l | 2 

= i <oi $| +) (x B ,t B )$W(xA,tA) m A )) l i 2 

= \A(x A ,t A ,x B ,t B )\ 2 (A4) 

where Ui is the interaction picture time-evolution oper- 
ator, the interaction picture field operator is 

$| +) (x, t) = V~ 1/2 «k,a exp[i(k ■ x - ckt)} (A5) 

k . a 

and in passing to the second line we have made the as- 
sumption that a detection at (xa, £a) does not physically 
influence the behavior of photons at (xs,ts) although 
there may be correlations. The the amplitude of joint 
detection is 

^(xA,iA,x B ,t B ) 

(27T) 3 



v 2 h 



Y^exp[i(k' • xa — ck'tA)] 



k'l' 



x exp[i(l' • x B — cl'tB)} 
x^ 5 (k)/i(l) W k 'i k 

kl 

1 - cxp[ic(fc' +V - k- l)t A ] 
c(k' + 1' - k-l) + irj 



(A6) 
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Although Vi'k'ik depends strongly upon the magni- 
tudes of the initial and final photon momenta through 
the resonance denominators of Eq. ([lO]) , it depends only 
weakly on their directions through the geometrical fac- 
tors of Eq. (^|). This and the assumption of paraxial 
input beams justify the approximation 

^2g{k)h(l)V VVSk 

kl 

~ Vvk>\ k ^2 g(k)h(l)5k + i t v + v 

kl 

= Vi'k'ioko / d 3 xG{x)H(x) exp[-i(k' + I') • x]. (A7) 

We can similarly treat the output photons in the parax- 
ial approximation for the case that the detection points 
are far from the interaction region, i.e., that xa, x b ^> x. 
Making these approximations and dropping unphysical 
portions of the solution propagating inward from the de- 
tectors toward the source region, we find 

A(-XA,tA,*-B,t B ) 

= J— I fe'dfc'/Vi'k'loko 



d A x- 



G(x)#(x) 



\x A - x||x B - x| 
x exp[z(k' • (xa — x) — ck'tA)] 
x exp[z(l' • (xb — x) — cl't B )} 
*0{t a )9(t b ) 



(A8) 



where cta.b = ctA,B ~ %a,b are retarded times. 

A final approximation ignores the slow variation of 
k', I' relative to that of the resonant Vi'k'i k - Further, we 
define G'(x) = G(x) exp[ik -x], H'(x) = H(x) exp[ilo -x] 
and k' = k' + 8k' where kp is the value of k' which max- 
imizes Vi'k'i k subject to momentum and energy conser- 
vation. This gives a simple expression for the correlation 
function 

exp[-ic(k' T A + 1' t B )] 



Tic 



x J d5 k ,V V ui K exp[ic<5 fc , (r B - t a ) 
G\x)H'{x) 



x / <Fx- 



\x A - x\\x B - X| 
x exp[i(k + 1 - ko - 1' ) • x] 
x0(t a )6{t b ). 



(A9) 



This can be interpreted as consisting of a carrier wave, 
a Fourier transform of the scattering amplitude and a 
coherent integration of the contributions from different 
parts of the interaction region. The spatial integral en- 
forces phase matching in the photon-photon scattering 
process. 



2. Detection Rates 

The probability for a coincidence detection at two de- 
tectors of specified area and in two specified time inter- 
vals is 

P = J d 2 XAd 2 x B cdtAcdt B \A(xA 1 t A ,x B ,t B )\ 2 , (A10) 

where the integral is over the detector surfaces (each as- 
sumed normal to the line from scattering region to de- 
tector) and over the relevant time intervals. This is more 
conveniently expressed in terms of a rate W of coinci- 
dence detections in terms of the detector solid angles 
S£Ia, Sfl B and the difference in retarded arrival times 

T- = T B - T A 



w 



c x 



v x 2 B \A(TC A ,t A ,x B ,t B )\ 2 6n A SQ B dT-. (All) 



Coincidence rate is largest when the detectors are 
placed in the directions which satisfy the phase-matching 
condition. We assume that k + l = k' + l' = and that 
the detectors are small compared to the source-detector 
distance, i.e., that SHa.e <C 1. Under these conditions, 
the rate of coincidence events reduces to 



scattering 



(k'l 



/7/\2 



J dSk'V\'k'i k exp[icJfc'T_ 



d 3 xG(x)H(x) 



5Q, B 8VL B dT-. (A12) 



3. Signal Contrast 

In addition to the photon-photon scattering signal, un- 
corrected Rayleigh and Raman scattering events give a 
background of accidental coincidences. The rate of scat- 
tering into a small solid angle 5Q is 



w BG = ssn 



xn k 



(A13) 



where 



B = J2 I d 3 p(n a , p )(l±(n c , p ,)) — 

a,c 

><E 



kjc 



(2-Kfti 2 

0/ • Vbc)* e i ■ Mba 



ck - 



Jab- ^[p-k + fc 2 /2] +ij b 



(A14) 



and rife is the number density of photons propagating in 
the k direction. In terms of the beam-shape functions for 
two colliding beams, the rate of accidental coincidences 
is 



G 



^accidental = B 2 j d 3 x|G(x)| 2 + \H(? 
xSCIaSCIbcIT-. 



(A15) 



The ratio of coincidences due to photon-photon scatter- 
ing to accidental background coincidences is thus 



^scattering _ jk'l') 2 F 
^accidental 4ft 2 B 2 

x J d5 k > Vi'k'i k cxp[ic4/r_ 

where F is the mode fidelity factor 
[fd 3 zG(x)tf(x)] 2 



F = 4- 



[Sd?x{\G{*)\ 2 + \H{*)\ 2 )Y 



(A16) 



(A17) 
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FIG. 1. FIG. 1. Photon-photon scattering at a single atom. 




FIG. 2. FIG. 2. Right-angle scattering amplitude A vs. time delay for the single-atom process of Fig. 1. The time unit is 
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FIG. 3. FIG. 3. Right-angle scattering rate |^4| 2 vs. time delay for the single-atom process of Fig. 1. Time unit is S^ 1 . 




FIG. 4. FIG. 4. Statistical correction to the one-atom scattering amplitude. The time unit is S^. 
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FIG. 6. FIG. 6. Scattering rate \A\ 2 vs. time delay for the two-atom process of Fig. 5. The time unit is 5^. 
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FIG. 7. FIG. 7. Scattering rate \A\ 2 vs. time delay for a Bose gas of phase-space density 1/2. The time unit is S^, 1 . 



8 




-4 -2 2 4 

TB-^A 

FIG. 8. FIG. 8. Scattering rate \A\ 2 vs. time delay for a Fermi gas of phase-space density 1/2. The time unit is S^, 1 . 
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FIG. 9. FIG. 9. "Ladder" process in a three-level atom. 




FIG. 10. FIG. 10. Scattering rate \A\ 2 vs. time delay for the "ladder" process of Fig. 9. The time unit is 8^. 
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FIG. 11. FIG. 11. Scattering rate \A\ 2 vs. time delay for a Lorentz-model atomic medium. The time unit is S^ 1 . 
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FIG. 12. FIG. 12. Coincidence rate vs. time delay for a Lorentz-model atomic medium. The constant background is 
accidental coincidences due to independent Rayleigh scattering events. The time unit is S^, 1 . 
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FIG. 13. FIG. 13. Small angle scattering rate \A\ 2 vs. time delay for nearly co-propagating photons. Time unit is S^, 1 . 
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FIG. 14. FIG. 14. Right-angle scattering rate \A\ 2 vs. time delay for nearly counter-propagating photons. Time unit is <5 ; 
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FIG. 15. FIG. 15. Coincidence rate \A\ 2 vs. time delay for nearly co-propagating photons in a linear medium. j b — O.Olfo, 
A = 2-kS-d. Time unit is S^ 1 . 
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FIG. 16. FIG. 16. Coincidence rate \ A\ 2 vs. time delay for nearly counter-propagating photons linear medium, = O.Olfo, 
A = 2-kS-d- Time unit is S^ 1 . 



10 




11 



